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ABSTRACT 

We combine new VLA D array HI data of NGC 2613 with previous high 
resolution data to show new disk-halo features in this galaxy. The global HI 
distribution is modeled in detail using a technique which can disentangle the 
effects of inclination from scale height and can also solve for the average volume 
density distribution in and perpendicular to the disk. The model shows that 
the galaxy's inclination is on the low end of the range given by Chaves & Irwin 
(2001) and that the HI disk is thin (zg = 188 pc), showing no evidence for halo. 

Numerous discrete disk-halo features are observed, however, achieving z 
heights up to 28 kpc from mid-plane. One prominent feature in particular, of 
mass, 8 x lO'^ M0 and height, 22 kpc, is seen on the advancing side of the galaxy 
at a projected galactocentric radius of 15.5 kpc. If this feature achieves such 
high latitudes because of events in the disk alone, then input energies of order 
~ 10^^ ergs are required. We have instead investigated the feasibility of such a 
large feature being produced via buoyancy (with drag) within a hot, pre-existing 
X-ray corona. Reasonable plume densities, temperatures, stall height (~ 11 kpc), 
outflow velocities and ages can indeed be achieved in this way. The advantage of 
this scenario is that the input energy need only be sufficient to produce blow-out, 
a condition which requires a reduction of three orders of magnitude in energy. If 
this is correct, there should be an observable X-ray halo around NGC 2613. 

Subject headings: galaxies: individual (NGC 2613), ISM: bubbles, galaxies: halos 



1. Introduction 



High latitude HI in disk galaxies is increasingly being recognized as an important tracer 
of disk activity and dynamics. The gas scale height is a measure of time integrated activity 



- 2 - 



in the disk and discrete features trace the level of local activity. The disk-halo interface 
represents a severe transition between two very different environments: the high density, low 
velocity dispersion disk and the low density, high velocity dispersion halo. Globally, it may 
provide information on the shape of the dark matter halo (e.g. Oiling & Merrifield 2000). 
In a few recent cases, lagging HI halos have also been detected, such that the high latitude 
gas is rotating ~ 20 - 50 km s~^ more slowly than the underlying disk (Schaap, Sancisi, 
& Swaters 2000, Praternah et al. 2001, Lee et al. 2001, Rand 2000 and Tiillmann et al. 
2000). Explanations for such lags have focussed mainly on galactic fountains (cf. Bregman 
1980, Swaters et al. 1997, Schaap et al. 2000, CoUins et al. 2002). Edge-on galaxies are 
particularly useful in the study of the disk-halo region and, in this paper, we especially wish 
to investigate how discrete disk-halo features can be used as probes of the environment into 
which they are emerging. 

The edge-on galaxy, NGC2613 (D = 25.9 Mpc) is a good candidate for such a study. In 
Paper I (Chaves & Irwin 2001), we presented Very Large Array (VLA) C Array observations 
of this galaxy. The resulting C array zeroth moment map over an optical image is shown 
in Fig. 1, illustrating the presence of six symmetrically placed vertical extensions in this 
galaxy. In this paper, we present new VLA D array HI data of NGC 2613 and examine 
them separately as well as in combination with the previously obtained C array data. The 
observations are presented in Sect. 2, the results, including a description of extremely large 
high latitude HI features and a model of the global HI distribution are given in Sect. 3, in 
Sect. 4, we present the discussion, including a feasibility study of a buoyancy model, and a 
summary is given in Sect. 5. 



2. Observations 

Details of the C array data acquisition and mapping have been given in Paper I. 

The D array data were taken on Sept. 15, 2000. The flux calibrator was 3C 48 and the 
phase calibrator was 0837-198. The flux calibrator was also used as a bandpass calibrator and 
was observed once for 17 minutes at the beginning of the observing period. The phase cali- 
brator was observed in 3 minute scans alternating with 25 minute scans on NGC 2613. The 
total on-source observing time was 2.4 hours and on-line Hanning velocity smoothing was 
applied. The data were edited, calibrated, and Pourier Transformed using natural weighting, 
cleaned, and continuum subtracted in the usual way using standard programs in the Astro- 
nomical Image Processing System (AIPS) of the National Radio Astronomy Observatory. 
The resulting (RA, DEC, Velocity) cube was then examined for broad scale structure. 
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The calibrated D array UV data were then combined with the previously obtained C 
array UV data (on-source observing time = 15.3 hours, see Paper I) to improve the signal to 
noise (S/N) ratio. Since the high resolution C array data have already been presented, we 
here concentrate on the broader scale structure by re-mapping the C+D array data, using 
natural weighting, over a UV range corresponding to the D array data only, i.e. 5 kA. That 
is, we show the lower resolution results (at a typical D array resolution), but with the highest 
S/N possible. Table 1 lists observing and map parameters for D array and for the combined 
C+D array data. Note that, although we show images for the combined data set, we have 
inspected and analyzed all data sets both independently and in combination. 



3. Results 

Note that a full description of the C array results and relevant images have been given 
in Paper I and we have reproduced the zeroth moment C array map in Fig. 1. The remaining 
The zeroth and first moment maps from the combined data set are shown in Fig. 2, both 
rotated so that the x axis is parallel to the major axis of the galaxy. Vertical shces have 
been taken perpendicular to the major axis at positions numbered 1 through 10 in Fig. 2a 
and these are shown as position-velocity (PV) plots in Fig. 3. The PV slices in the last two 
panels of Fig. 3 show averages over the receding and advancing sides of the galaxy. 



3.1. Tidal Tail 

For the first time, we see that there is a tail of emission on the south-eastern (left. 
Fig. 2a) side of the galaxy which resembles tidal tails seen in galaxy interaction simulations. 
This tail is likely produced via an interaction with the companion, ESO 495-G017, to the 
north-west, whose systemic velocity is separated from that of NGC 2613 by only 143 km s~^. 
The tail is at negative velocities with respect to the receding side of the galaxy from which 
it originates and shows a strong velocity gradient along its length, becoming negative with 
respect to the galaxy's systemic velocity as well, i.e. the velocities become forbidden. This 
can occur in the case of interactions but is difficult to explain any other way. Given this 
velocity structure, the tail must be traihng from the eastern, receding side of NGC 2613 and 
is therefore in front of the galaxy. The anomalous emission seen around -200 km s~^ in Slice 
2, as well as the average over the receding side of the galaxy (second last panel) of Fig. 3 
belongs to this tail as well. There is excess emission on the north-west (right. Fig. 2a) side 
of the galaxy in the vicinity of the companion, but wc have foimd no systematic velocities 
associated with that emission or any other sign of a corresponding redshifted feature. 
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3.2. High Latitude HI 

In Paper I, using the C array data alone, we identified and characterized the properties 
of 6, symmetrically placed extensions labelled Fl to F6 in Fig. 1. Fig. 2a, which includes 
the lower D array resolution data, now also reveal these same features on broader scales, 
some of which are blended with other features. Since a single cutoff level must be applied 
over the whole galaxy to create this image, this moment map is best used to see a global 
view of the disk-halo features and identify where new features might be present. The full z 
extents and, of course, the velocities of the individual features are best measured from PV 
slices perpendicular to the major axis. We have inspected many PV slices over the map 
and show selected slices (labelled 1 through 10 in Fig. 2a) in Fig. 3. In the new data, the 
previously seen features Fl through F6 occur at the same velocities and the same positions 
as shown in Fig. 1, after taking into account the different beam sizes and the addition of 
broader scale emission. Please refer to both Figs. 2a and 3 for a more thorough discussion 
of these features, presented below. 

Fl appears similar to what has previously been seen (Fig. 1), but at lower resolution. 
The PV plot for Fl can be seen in Slice 2 which reveals it as a distinct feature at a velocity 
of 250 km s~^ but it can also be seen in Shce 1 as well, extending to a projected z height 
of at least 100" (12.6 kpc) "above" (i.e. positive z) mid-plane. A newly seen extension 
immediately to the west of Fl, sampled by Slice 3, may actually be associated with or part 
of Fl, since it is seen in PV space as a small arc at about the same velocity as Fl. 

F2 (see Slice 1) shows several velocity components and is more complex than previously 
seen. Like Fl, F2 extends to at least 100" at a velocity of 210 km s~^. It is not clear whether 
the other components may be blended with emission from the tidal tail. 

F3 reveals more structure than the extension in Fig. 1 and is sampled by Slice 4 which 
reveals it as a 2- pronged arc above the disk centered at a velocity of 100 km s~^. Such double 
peaks in velocity space are typical of such features and have been seen in other galaxies (see, 
e.g. Lee & Irwin 1997, Lee et al. 2001). 

The extension to the left of the F4 label in Fig. 2a) extends to a height 225" (28 kpc) 
below (i.e. negative z) mid-plane in this map. As can be seen in Slice 3 of the PV map, this 
feature occurs at a velocity of 50 km s~^ and shows no velocity gradient with height. F4 
reaches the highest ^of all the extensions and, to our knowledge, reaches higher latitudes than 
any previously seen in an edge-on galaxy, exceeding even the starburst galaxy, NGC 5775, 
whose known HI features extend to ~ 7 kpc from the plane (Lee et al. 2001). This is 
remarkable, given that NGC 2613 (log \jpir = 10.00) has a massive (0.1 ^ 100 Mq) star 
formation rate that is a factor of 3 lower than that of NGC 5775 and a supernova input 
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energy rate per unit area that is a factor of 7 less than NGC 5775 (Irwin et al. 1999). 

Feature F5 is sampled by Slices 8 and 9 and can most clearly be seen in Slice 9 as 
a vertical feature above the plane in PV space at a velocity of -300 km s~^. Whereas it 
appeared as a single extension in Fig. 1, it can now be seen to have a counterpart 55" (7 
kpc) in projection farther out along the disk which is sampled by Slice 10. The counterpart 
clearly occurs at the same velocity and also shows no velocity gradient with height. These 
two features may form the base of a large loop or cylinder, similar to what has been seen in 
galaxies hke NGC 5775 (Lee et al. 2001). 

The most interesting and well-defined feature, however, is associated with F6. For 
example, slices 6, 7, 8 and 9 pass through a feature below the major axis which extends 
22 kpc (175") in projection, from the plane and centers at -307 km s~^. Individual channel 
maps of this feature are shown in Fig. 4. We will refer to this as the -307 km s~^ feature, 
rather than F6 since, as revealed by the new data, the feature is more extensive than what 
was originally seen in Fig. 1. This feature will be discussed extensively in Sect. 3.4. 

This disk-halo emission is faint (3 - 4 o" typically, note contour spacings in Fig. 3) but 
we consider the above-discussed features to be real because, 1 ) their spatial morphology is 
typical of the kinds of extensions and loops seen in other edge-on systems, 2) even though the 
feature might be faint far from the disk, there is evidence for disturbed contours 'beneath 
them' even at highly significant contour levels (see the 5th and 6th contours of Fig. 2a, 
for example), 3) at least the 6 major features, Fl to F6 can be seen in both independent 
data sets, 4) the features generally occur over several independent beams and/or velocity 
channels, and 5) in PV space, the features arc distinct, yet smooth and connect continuously 
to the disk, unlike artifacts produced by a badly cleaned beam (see, e.g. the arc above the 
disk seen in Slice 3 of Fig. 3 at ~ 270 km s^^ or the larger structures in Slices 9 and 10). In 
addition, 6) the Receding and Advancing panels of panels of Fig. 3) show averages over the 
entire advancing and receding sides of the disk, rather than just selected individual slices as 
shown elsewhere in the figure. In these panels, therefore, we have basically averaged together 
signal with noise and weak extensions would be diluted (i.e. "beam-averaged", where the 
beam is the size of half the galaxy) . Yet these panels still show the extensions quite clearly, 
especially on the Advancing side where the disk-halo emission is dominated almost entirely 
by the -307 km s~^ feature. This feature is seen up to the 4th contour which is a 6(7 detection 
in the diluted averaged panel. Finally, 7) for at least 4 of the 6 features, there is evidence 
from independent observations that extraplanar emission is present. For example, our radio 
continuum image (see Fig. 9 of Paper I) shows extended emission at the positions of F3, F4, 
F5, and F6. Also, an early HI image by Bottema (1989) at lower resolution shows extended 
HI emission at roughly the same locations. 
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The zeroth moment map of Fig. 2a also shows a number of low intensity "disconnected" 
features at high latitudes (for example, the emission around -300", -200"). We do not claim 
that all of these are real, but we have chosen to display them for the sake of comparison with 
future observations and because some of the emission shown in independent lower resolution 
images (e.g. Bottema 1989, Irwin et al. 1999) is so extended. An example would be the 
'knot' of emission between F3 and F5 which is sampled by Slice 6. For this feature, the 
velocities are discordent with respect to the underlying disk (see Fig. 3) and therefore do 
not satisfy the criteria hsted above. 

Particularly interesting about these very high latitude features is that there is little or 
no evidence for any lag in velocity with z, even to very high distances from the place. The 
features seen below the plane in slices 3, 8, and 9 and features seen above the plane in slices 
9 and 10 show no velocity gradient. The early data of Bottema (1989) also showed that the 
high latitude HI was co-rotating. There is some evidence, in our data, for curvature in the 
smaller features (cf. Slices 3 & 4 above the plane), but one could not make any case for a 
lag. The only slice which shows some evidence for lagging gas is Slice 1, but this slice is 
also closest to the tidal tail, confusing the interpretation. The averages of the advancing and 
receding sides of the galaxy also do not support the presence of a lagging halo. Here, we see 
an asymmetry in the sense of many disk-halo features on the receding side, but the advancing 
side is dominated by the single -307 km s~^ feature. None of the features, however, show a 
decline with z. We investigate this further in the next section. 



3.3. Kinematic Models of the Global HI Distribution 

We have modeled the global HI density and velocity distributions for NGC 2613 following 
Irwin & Seaquist (1991) and Irwin (1994). This approach models all spectra in the HI cube 
by adopting a volume density distribution in and perpendicular to the plane, a rotation 
curve and velocity dispersion (if required) , and a position and orientation on the sky. Given 
trial input parameters describing these distributions, the routine then varies the parameters, 
examining the residuals, until the best fit solution is found. Various sections of the galaxy 
can be modeled, including galactocentric rings in the galaxy's plane. In the current version, 
we can also isolate slices parallel to the galaxy's plane. This permits modeling the halo, 
independent of the disk, which is important in this study in the event that a lagging HI halo 
might be present. While there is no evidence for lags in the discrete features shown in Fig. 3, 
it is possible that the smooth, low intensity halo might be better fit with a lagging velocity 
distribution, similar to what has been found for NGC 2403 (Fraternali et al. 2001). The 
routine, called GUBIT, interfaces to classic AIPS; potential users are requested to contact 
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the first author for the code. 

For NGC 2613, the adopted in-plane density distribution is a gaussian ring, described 
by: n(R,z = 0) = noexp((R — Rof/D"^), where R is the galactocentric radius, Rq is the 
galactocentric distance of the center of the ring, uq is the in-plane density at the center 
of the ring, and D is the scale length. D can have different values for points R > Rq 
{Do) and for points R < Rq (Di). The distribution perpendicular to the plane is given by 
the exponential, n{R,z) = n{R,z = 0) exp{~z/Hf,), where z is the perpendicular distance 
from midplane and is the exponential scale height. The rotation curve is described by 
the Brandt curve, V{R) = Vsys ± Vmax (R/Rmax) {'^/S + 2/3 (i?/i?„,aa;)'")^~^^^™\ where 
Vsys is the galaxy's systemic velocity, Vmax is the peak of the rotation curve occurring at 
a galactocentric radius of Rmax and m is the Brandt index which is a shape parameter. 
A gaussian velocity dispersion, ay can also be applied, if needed. Note that a velocity 
dispersion will result from all contributions to line widening along the entire line-of-sight 
direction, for example, any non-circular motions along the plane, streaming motions due to 
spiral arms and local turbulence. Thus, cry may be larger than what might be expected from 
a local value. The above distributions result in 10 free parameters: no, -Ro, Do, Di, He, Vgys, 
Vmax, Rmax, ""^j a^d (Ty. In addition, there are 4 more orientation parameters: the position 
of the galaxy center, RA(0), DEC(O), the major axis position angle, PA, and the inclination, 
i. Over the region of the cube occupied by NGC 2613, there are 552 independent data points 
(beam / velocity-channel resolution elements) above a 2 cr intensity level in the combined 
C-l-D cube, so the 14 parameters are well constrained, over all. Certain parameters may 
show a larger variation, however, depending on the nature of the distribution; for example, 
Rmax tends to have a larger error since the point at which the rotation curve changes from 
rising to flat is not as well constrained by the data as some of the other parameters. 

Table 2 hsts the best fit modeled values for NGC 2613 for the combined C-l-D array 
data. We also modeled the C array data and D array data alone; the error bars of Table 2 
reflect the variations in the parameters over the different arrays that result. Given the range 
of resolution over these arrays, from 26.1" x 20.2" (3.3 x 2.5 kpc) to 96.5" x 42.9" (12.1 
X 5.4 kpc) the error bars on the modeled parameters are very small. The best fit model is 
shown as dashed curves superimposed on the data in Fig. 3. Taking only residuals (data 
cube - best fit model cube) greater than 3a, the root- mean-square value is 1.77 mJy beam~^ 
and the average relative error is ~ 10%. A zeroth moment map made from the cube of the 
residuals of the best fit is shown in Fig. 5. 

To determine whether there are changes in i, PA, and with galactocentric radius, R, 
we also modeled the galaxy in concentric rings. For example, either a lower value of i with 
increasing R (a warp) or a higher value of He with increasing R (i.e. a fiare) could conspire 
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to make the halo of a galaxy appear to lag. We find that the inclination does not decrease 
with radius (there is a slight increase). There is indeed evidence for flaring such that, in 
each data set, Hf. systematically increases with R. The most significant change is measured 
for the D array data alone such that the mean scale height between a radius of 140" and 
300" (17.6 to 38 kpc) is 5.5" (693 pc) compared to a value of 1.7" (214 pc) over the whole 
galaxy. Thus, the disk is quite thin with a modest fiare at large galactocentric radii. 

In order to model the halo alone, a variety of lower and upper cutoffs to z were applied 
and these data were modeled in a similar fashion. (Note that z refers to the actual distance 
from the galaxy's midplane, not just the projected distance from the major axis.) The 
number of independent data points is reduced typically by about a factor of 3, depending on 
which z limits are chosen, but we have also held some of the parameters (the central position, 
position angle and systemic velocity) fixed. For the halo models, we find no evidence that 
Vmax could be lower at higher z. The same conclusion is reached if the C array and D array 
data are modeled independently. The absence of a lagging HI halo in NGC 2613 will be 
discussed in Sect. 4.1. 



3.4. The -307 km s"^ Feature 

A distinct feature occurring at -307 km s~^ (see the Advancing panel of Fig. 3 and 
Fig. 4) has been mentioned in Sect. 3.2 and extends to a (negative) z height of 22 kpc. 
(This feature is also distinguishable on the major axis PV slice of Paper I). The channel 
maps show a complex structure below the plane, although two spurs dominate. The feature 
is very well defined in PV space and contributions to the feature can be seen in slices 6 
through 9 which covers a range of 13 kpc in projected distance along the disk of NGC 2613 
(Fig. 2a). There is also a distinct height {z ~ 90" = 11 kpc) at which the feature begins to 
widen abruptly in velocity (Fig. 3) and also (approximately) spatially (Fig. 4), a point we 
return to in the model of Sect. 4.3. Further support for HI at high z at this position comes 
from early HI observations by Bottema (1989) which show an HI "blob" extending to ~ 25 
kpc from the plane and sharing the rotation of the underlying disk. In addition, our low 
resolution 20 cm continuum map also shows a single large feature at this position extending 
21 kpc from midplane (Irwin et al. 1999). Thus, we will treat the -307 km s~^ feature to be a 
single, coherent structure, as implied by its appearance in PV space, centered at a projected 
galactocentric distance of 15.5 kpc. 

This feature occurs over a restricted region of observed radial velocity (K = -307 ± 30 
km s~^). Assuming that the feature originates in the disk (see Sect. 4.2), this should allow us 
to place limits on the position along the line of sight from which it originates since the rotation 
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curve of the galaxy is known. The galactocentric radius is given by R — Rproj {V{R)/Vr) 
(ignoring sin i = 0.98) where Rproj is the projected galactocentric radius and V{R) is the 
circular velocity at R which is approximately constant in this part of the galaxy (Table 
2). Allowing a ±30 km s^^ range on Vr, the feature could therefore have originated from 
R = Rproj = 15.5 kpc or anywhere within ± 7 kpc from this value along the line of sight, 
corresponding to a range of R between 15.5 and 17 kpc. Since the projected size of the 
feature is larger than this range (Fig. 4), it is likely that there has been some expansion with 
distance from the plane. 

Wc can also put a limit on the age of the feature. Gas which leaves the disk at an earlier 
time shouki have a component of line of sight velocity corresponding to the motion of the 
underlying disk at that time. As the disk rotates, gas which leaves later may have a different 
line of sight velocity component. Consequently, there may be a gradual change in observed 
velocity with height (increasing or decreasing) depending on geometry. Since we see no such 
change, we can use the velocity half- width of the feature, 30 km s~^, to place a limit on the 
angle over which the underlying disk has rotated over the age of the feature. We find that the 
feature must have formed over a timescale which is less than the time required for the disk 
to rotate ± 24° at galactocentric radii between 15.5 and 17 kpc. From the circular velocity 
in this region, the corresponding timescale is 4.2 x 10^ yr. These arguments will depend 
on how the ejected gas may or may not be coupled to the disk and should be considered 
order of magnitude only. Nevertheless, the result is similar to the kinematic ages found for 
expanding supershells in other galaxies (cf. King & Irwin 1997, Lee & Irwin 1997, Lee et al. 
2001). Basically the "straightness" of the feature in PV space requires that the feature have 
an age << the rotation period. 

The mass of this feature, as measured in PV space to the level at which the feature 
blends with the disk, is (8 ± 2) x 10^ Mq. Approximating the feature as a cylinder of 
height, 22 kpc, width and depth, 13 kpc, the mean density is 1.1 x 10~^ cm~^. (Reducing 
the diameter of the cylinder to 7.3 kpc (150") results in a mean density of 3.6 x 10~^ cm~^.) 
If the feature has reached a z height of 22 kpc in 4.2 x 10^ yr, the implied mean outflow 
velocity is 512 km s~^. The corresponding HI mass outflow rate would be 1.9 Mq yr~^. 

Wc can also compute the potential energy of this feature (see Lee et al. 2001) which 
requires a knowledge of the mid-plane stellar density and scale height. Comparing the 
rotation curve of NGC 2613 (Chaves & Irwin 2001) to that of the Milky Way (Sofue & 
Rubin 2001), NGC 2613 is a factor of 4 more massive, and factor of 2 larger in radius. 
NGC 2613 also has a thinner disk, if we take the thickness of the HI disk (188 pc. Table 2) 
to be representative of the stellar thickness. Scahng the mid-plane stellar density from the 
Galactic value (0.185 M© pc-^) yields 0.74 M© pc-^ for NGC 2613 and adopting 100" (12.6 
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kpc) (about half the full z extent) to be representative of the height of the feature above the 
plane, we find a potential energy of 1.5 x 10^^ ergs for the -307 km s~^ feature. Thus, energy 
of this order is required to transport the cool gas to the heights observed. Although very 
large, this value is not unlike that determined for the largest supershells in other edge-on 
galaxies (cf. an input energy of ~ 3 x 10^^ ergs for NGC 3556, King & Irwin 1997). We 
return to this issue in Sect. 4.3. 



4. Discussion 

4.1. The Absence of a Lagging HI Halo in NGC 2613 

There are now several galaxies for which there is evidence for lagging halos. In NGC 891 
(i > 88.6°) , the rotation curve of the HI halo between ^ = 1 .4 to 2.8 kpc reveals a peak velocity 
~ 25 km s-^ lower than in the disk (Swaters et al. 1997). In NGC 5775 (i = 86°), the ionized 
gas velocities decrease with z to heights of 5 - 6 kpc above which the velocity remains roughly 
constant (Rand 2000). The HI in this galaxy shows more complex structure, though lagging 
gas tends to dominate at z heights up to 5 to 6 kpc. Above this region, the velocities are 
more nearly constant or the features break up into clumps which are seen over a wide range 
in velocity (Lcc ct al. 2001). Models of NGC 2403 (inclination of 61°) suggest that a lagging 
halo extends to z ~ 3 kpc (Schaap et al. 2000). Thus, from the sparse data available, the 
region of the lag extends to typically 3 to 6 kpc above the plane above which (if gas is 
detected at all) there isn't strong evidence for a global lag. Given the presence of lagging 
halos in these other galaxies and also some theoretical expectation of lags (cf. Bregman 
1980, Benjamin 2002, CoUins et al. 2002), the absence of a lagging velocity gradient along 
the discrete features (Sect. 3.2) and the absence of a global lag (Sect. 3.3) in NGC 2613 
require comment. 

Firstly, given the inclination of this galaxy (79°, Table 2) and its HI radius of ~ 35 kpc, 
the discrete features must reach a height of z > 6.5 kpc before they will be seen beyond 
the projection of the disk. If it is true that lags only exist up to 3 to 6 kpc and above 
this there is no longer a lag (see above), then lags (if they exist) along the discrete features 
would be projected against the background disk. This would make a velocity gradient (if 
present) impossible to detect along the discrete features. The 'straightness' of the -307 km 

feature in velocity (Fig. 3) to such high latitudes beyond the projection of the disk is, 
however, quite remarkable. If this feature emerges into a hot X-ray corona, as suggested in 
Sect. 4.3, then either the velocity of the broad-scale X-ray corona is not appreciably different 
from the underlying disk (although the velocity dispersion is likely much higher), or the 
timescale over which it could appreciably affect this gas column is greater than the age of 
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the feature. 

As for a globally lagging halo, our kinematical model is indeed capable of detecting such 
a lag, even if it is projected against the disk (see Sect. 3.3) provided there is sufficient high 
latitude emission that halo gas can be detected at all. What we have found, however, is that 
NGC 2613 does not have an HI thick disk or halo. The vertical exponential scale height (188 
pc. Table 2) indicates that the disk is, in fact, thin. If we create a model galaxy using the 
parameters of Table 2, it is straightforward to show that emission from all high latitude (z 
> 1 kpc) HI in fact falls below the map noise. Thus, NGC 2613 does not have a detectable 
lagging halo because it does not have a halo at all. 

As a further comparative test, we consider the HI halo of another galaxy at the same 
distance, NGC 5775 (i = 86°, D = 25 Mpc), for which we detected an HI halo with exponential 
scale height of 9.14" (1.1 kpc) using the same technique (Irwin 1994). If we create a model for 
the lower inclination (i = 79°) NGC 2613 in which the vertical density distribution declines 
exponentially from its in-disk value with the same scale height as NGC 5775, we find that 
it could easily have been detected, even at the lower inclination of NGC 2613. 

This result illustrates the importance of using a model, rather than PV slices alone, in 
drawing conclusions about the possible existence of lagging halos. For example, the model 
can disentangle inclination effects from the effects of a thick disk. For NGC 2613, since the 
modeled inclination is now known to be on the low end of the range quoted in Chaves & 
Irwin (2000), our earlier conclusion of a lag with zcan now be largely explained by projection 
against the background disk. This is succintly illustrated via the dashed curves in Fig. 3 
which show that our thin disk model provides an almost perfect fit to the data. 

4.2. The Discrete Features: Internal Origin and Energies 

We have argued elsewhere (Paper I), largely on the basis of above/below plane sym- 
metry (Fig. 1), that the observed features are internally generated and represent outfiows. 
While these new data show that the galaxy is indeed interacting, it is unlikely that most 
of the observed features are produced via cloud impacts since the cloud would have to pass 
completely through the disk, forming similar structures on both sides, in contrast to what 
is expected theoretically (see Santillan et al. 1999). Anomalous velocities are also more 
likely in the case of impacts, whereas the velocities observed in NGC 2613 are typical of the 
underlying disk. Thus, the new data are consistent with the outflow interpretation. We can- 
not, however, rule out the possibility of some impacting clouds in this system and it is also 
possible that the interaction may assist in the process of disk-halo dynamics, for example 
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via stimulating a starburst or instabilities such as the Parker instability. 

If the features are internally generated, their very large sizes {z up to 28 kpc, Sect. 3.2 
imply exceptional energies (~ 10^^ ergs) which are difficult to conceive of, especially at large 
distances from the nuclear vicinity and in a galaxy which is not a starburst. This energy 
problem has been known since the first detections of Heiles Shells in the Milky Way (Heiles 

1979, 1984) and has generated suggestions as far ranging as multiple supcrnovae and stellar 
winds (Heckman 2001), gamma ray bursters (Efremov et al. 1998, Locb & Perna 1998), and 
jet bubbles (Gopal- Krishna & Irwin 2000) to help explain the high energies. In the following 
section, we suggest that the exceptional heights achieved by the features in NGC 2613 may 
be as much a result of the environment into which the features emerge, as to the input source 
itself. 



4.3. Buoyant Outflow 

The -307 km s^^ feature (Fig. 3, Advancing) is remarkable in its well defined structure, 
in the very high latitude (22 kpc) that it achieves, and in its sudden and dramatic widening 
and break up in velocity at a height of ~ 11 kpc. In individual slices (sec slices 8 and 9, 
in particular), this velocity widening gives the feature a mushroom- like appearance and is 
reminiscent of buoyant gas rising through a higher density medium. A precedent for this 
kind of behaviour is the 350 pc Galactic Mushroom discovered in Canadian Galactic Plane 
Survey data and which has been interpreted in terms of buoyant outflow (English et al. 
2000). 

Given this similarity, we here consider whether buoyant outflow could explain the ob- 
served high latitude discrete HI features in NGG 2613. We present the following as a fea- 
sibility study only to see whether the results provide a reasonable match to the data using 
realistic parameters. We consider a match to the -307 km s^^ feature only, at this time, since 
it is the clearest case amongst the disk halo features in NGC 2613; however, if the model is 
correct, it must clearly apply to the other features as well. A similar development has been 
presented by Avillez & Mac Low (2001) for smaller features like the Galactic Mushroom, 
but we here consider much larger scales (many kpc), and also include the effects of drag. 

The scenario envisioned is one in which a hot. X-ray emitting corona already exists 
around the galaxy, possibly set up via venting through previous fountain or chimney activity. 
Some event or events occur within the disk which are sufficiently energetic to produce blow- 
out. The HI is already entrained or swept up in some way by the time the outflow emerges 
into the halo. An initial velocity (i.e. the velocity at blow-out) could be present but is not 
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included here. Thus, the velocities achieved in the plume are a result only of buoyant forces. 
The attraction of this model is that, rather than requiring that all of the energy in the 
HI plume be supplied by the instigating event in the disk (10^^ ergs. Sect. 3.4), we require 
only enough energy to produce the blow-out condition (e.g. 10^^ ergs, Tomisaka 1998). The 
remaining energy is extracted from the pressure gradient in the hot coronal gas. Ultimately, 
the galaxy's potential itself is providing the energy source. 

Although we do not investigate the details of the entrainment or sweeping up of HI, we 
assume that the hot gas inside the plume carries the cool HI with it and that the density 
of the HI declines with z in a fashion similar to the hot buoyant plume material. In this 
feasibility study, we do not consider structure in the corona or outflow plume, other than 
the cylindrical geometry chosen for the plume, and also neglect the effects of shocks. The 
integrity of the HI in the presence of hot gas will be considered in the next section. 

The stellar/mass density distribution of the thin disk is described by: 

p*{z) = p*(0) sec/i^(— ) 

where is the scale height of the stellar disk, and p*(0) is the stellar density at mid-plane. 
Integration of Poisson's Equation: 

^ = -47rGp^{z) 
az 

over the regime, z » z^, together with the above density distribution, results in a gravi- 
tational acceleration which is constant with z and has a magnitude: 

g = 8TiGp,{Q)z, (1) 

The hot coronal gas (subscripted, c) is taken to be isothermal at temperature, T^., with an 
exponential fall-off in both pressure and density, respectively: 

Pc{z) = Pc{z.)e-'"' p,{z) = pc{z.)e-'"' (2) 

where H is the scale height. This hot coronal distribution starts at the top of the (HI -|- 
stellar) thin disk, i.e. at z* where the outflow just blows out of the thin layer. The scale 
height is given by 

H = kTj{gpmp) (3) 

where k is Boltzmann's constant, p is the mean molecular weight and nip is the mass of the 
proton. We consider a pure hydrogen gas {p — 1) for simplicity. 
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Inside the plume (subscripted, i), we consider the gas to be adiabatic. The temperature 
gradient is: 

dTi , 1 Ti dPi 

— ^ = (1 ) — — - (4) 

dz ^ ^i' P, dz ^ ' 

For a plume in pressure equilibrium with its surroundings, 

Pi{z) = Pe(^*) e-^/^ (5) 

where, again, we take the plume base to occur at z^. Substituting Eqn. (5) and its derivative 
into Eqn. (4) yields: 

T,{z) = Ti{z.)e^^^^ (6) 
and from Eqn. (5), Eqn. (6) and the perfect gas law: 



Pi{z) = pi{z,)e (7) 



The equation of motion (force per unit volume) of the plume material is: 



Pi{z) ^ = pc{z) g - pi{z)g - CdPc{z) ^ - (8) 



The first term on the right hand side denotes the buoyancy force, the second term gives 
the weight of the plume material, the third term specifies the drag force against the upper 
surface (assuming a case in which the effects of this force can propagate through the column) 
where C/j is the dimensionless drag coefficient, and the fourth term represents the viscous 
drag on the cylinder sides, where r is the cylinder radius and 77 is the coefficient of viscosity 

Both Cd and 77 depend on the Reynolds number: 

i?e fa 5 X 10 



14 pc{z) LvlnA 



rp5/2 

assuming a pure ionized hydrogen gas, where L is a scale length above which motion is 
damped by viscous effects and 



A Pi 1.3 X 10^ 



j,3/2 



where ne{z) is the coronal electron density. Taking v ~ 500 km s"-*- (Sect. 3.4) and Tc = 2 x 
10^ K, we find for size scales, L, of order several kpc. Re is in the range 10^ In this range 
of Re, Cd ~ 1 for an incompressible fiuid with cylindrical geometry. The viscosity coefficient 
includes both turbulent and molecular terms, i.e. 77 = Vt + Vm, where 77^ = Lv Pc{z)/Re 
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and Tjt — {Re/ Recrit) Vm, where Rccrit is the critical Reynolds number which designates the 
value of Re at which the flow becomes turbulent. Re/Recru is not known and depends on the 
geometry and nature of the interface, but typically has values between 1 and 100. We can 
consider whether the drag force on the cylinder sides (4th term of Eqn. (8)) is appreciable 
in comparison to the drag on the top of the cylinder (3rd term of Eqn. (8)). Taking L ~ r, 
the ratio of the 3rd to 4th drag terms becomes [(-Re/2) (r/z) 1/(1 + Re/Rccrit)]- Since r is 
of order z, then for Re ~ 10^"*^, even Re/Rccru up to 100 ensures that the drag at the top 
will dominate over that at the sides. Therefore the 4th term in Eqn. (8) is small and will be 
neglected. 

Substituting v = dz/dt, 7l{z^) = Pc{z*) / Pi{z^) , and 7 = 5/3 into Eqn. (8) and rear- 
ranging yields: 



^ = n{z.) e-°-4i n{z.) e-o-^i - (9) 

dz V V z 

This equation was integrated numerically for the input parameters shown in Table 3, 
yielding the curves of v{z) and z{t) shown in Fig. 6 and Fig. 7, respectively. The stellar scale 
height, z^ is set to 188 pc (Table 2) and we consider a coronal temperature, Tc — 2.5 x 10^ 
K (Models 1 to 4) which is comparable to that found from X-ray observations of NGC 253 
(Pietsch et al. 2000) as well as a temperature which is a factor of 2 higher (Model 5). Since 
the midplane stellar density of 0.74 M© pc~^ is a value which has been scaled from Galactic 
values, given the size and mass of NGC 2613 (see Sect. 3.4), we also consider a slightly lower 
value (Models 3 and 4). The behaviour of the curves depends only on the density ratio at 
the base of the corona, TZ{z<^), rather than the individual densities, but we can additionally 
fix the density at the bottom of the corona, nc(-2*) to be equal to the HI density at the top 
of the thin disk at the position of the plume (see Table 2), providing constraints upon the 
density within and outside of the plume. The peaks of the velocity curves (Fig. 6) indicate 
where the acceleration of the phimc material goes to zero; we define the z height at which 
this occurs to be Zgtaii- This position should correspond to the point at which the plume 
widens in velocity space, observationally determined to be ~ 11 kpc. 

Models 1 and 2 (Fig. 6) show the effect of changing the initial density ratio, such that 
the lower ratio (Model 2) results in a lower maximum velocity and a lower Zstaii (7-4 kpc as 
compared to 11.9 kpc). A comparison of Models 1 and 3 or of Models 2 and 4 show the effect 
of decreasing the midplane mass density. This lowers the gravitational acceleration (Eqn. (1)) 
which lowers both the buoyancy and the weight of the plume material. It also increases the 
scale height (Eqn. (3)) which increases both the buoyancy and the drag. The net effect is 
higher Zgtaii at lower g for the range of parameters given here. A comparison of Models 1 
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and 5 shows the effect of increasing the coronal gas temperature. The higher temperature 
increases the scale height alone (Eqn. (3)) which, again, increases both buoyancy and drag. 
Since drag is velocity dependent, the net effect is that the peak velocities achieved are lower 
for higher coronal temperatures. 

Clearly, Models 1, 4, and 5 provide an adequate match to the observed Zgtaii- Model 1, 
however, results in a mean plume density (between z^, and Zstaii) of only 0.27 x 10~^ cm~^ 
whereas the HI density alone is ~ 1 x 10~^ (Sect. 3.4), and is therefore not realistic. Models 

4 and 5 both result in reasonable fits to the known observational parameters, with Model 

5 slightly preferred because of its higher internal plume density (higher than that of the HI 
alone) and lower internal temperature. Note that the mean internal temperature derived 
here (~ 10^ K) is comparable to the hot (T = 1.2 keV = 1.4 x 10^ K) component of the 
outflow in NGC 253 (Pietsch et al. 2000). The peak velocities derived here (~ 250 to 300 
km s~^) are lower than the gross estimate of 500 km s~^ computed in Sect. 3.4 since we 
only model the feature up to the stall height (11 kpc) rather than over its total length. The 
timescales of ~ 5 to 6 x 10^ yr (Fig. 7) are in good agreement with the estimate found from 
the kinematical structure of the plume (4 x 10^ yrs). The mass outflow rate to the stall point 
is ~ 1 Mq yr~^. Mass flow continuity requires that the plume of Model 5 should increase in 
radius by a factor of 1.8 between z heights of 4.7 and 11.1 kpc. 



4.4. Integrity of the HI 

It is unlikely that the HI will be uniformly distributed at its mean density of ~ 10^'^ 
cm^^ but may exist in denser clouds or clumps whose sizes and distribution are not described 
by this simple model. In general, however, the situation will be not unlike that of those 
high velocity clouds (HVCs) which are in the hot halo of the Milky Way and the same issues 
regarding whether or not the clouds can remain neutral must be considered. Photoionization 
of the HI by starhght should be neghgible, considering the high galactic latitudes achieved. 
As pointed out by Murali (2000), the relevant interactions are HI - proton and HI - electron 
interactions for which the interaction cross-sections are in the range a — 10~^^ to 10~^^ 
cm^^ for a relative velocity of order 200 km s~^ as is appropriate here. The mean free path 
into the HI cloud is then 1 = l/(ncO') = 0.3 to 3 pc or smaller if the HI is clumped. This is 
considerably smaller than the size scale of the plume and therefore these ionizing interactions 
should be minor. 

The most important interaction will be heating due to thermal conductivity leading to 
the evaporation of clouds. The classical mass evaporation rate, applicable to the case in 
which the mean free path is small in comparison to cloud size, is given, for approximately 
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spherical clouds, by m — 2.75 x 10^ T^/^ Rpc (f) g s"^, where T is the temperature of the 
external medium, Rpc is 1/2 of the largest dimension of the cloud in parsecs, and is a 
parameter which measures the inhibition of heat flux due to the magnetic field and cloud 
geometry (Cowie & McKcc 1979, Cowie et al. 1981). We will assume = 1 (no inhibition) 
which maximizes the evaporation rate. Using T = 2.5 x 10^ K and R = 5 x 10"^ pc we find 
fn = 0.021 Mq yr^^ (to within factors of a few, given the different geometry). The timescale 
for complete evaporation of the 8 x 10^ Mq HI plume at a constant rate is then M/ m — A 
X 10^ yr. Since this is two orders of magnitude larger than the age of the plume, we expect 
that the plume will not evaporate over its lifetime. A caveat, however, is that since M/ m 
oc R^, if HI is distributed in many smaller clouds, then the evaporation time could approach 
the age of the plume. 



4.5. The X-ray Corona 

If the HI is indeed reaching such high values of z because of buoyancy, then HI disk- 
halo features probe the parameters of the halo gas, as suggested in Sect. 1. We have so 
far considered only whether buoyancy is feasible. There may be other dynamics at work as 
well (for example, an initial velocity at blow-out or magnetic fields) and the proposed hot 
corona is also unlikely to have the smooth distributions postulated here. Nevertheless, in 
the context of the model, it is interesting to predict the X-ray luminosity of the corona in 
the vicinity of the plume. 

For the case in which the ion and electron densities are equal, the index of refraction is 
1, the charge Z =1, the frequency is greater than the plasma frequency, and T > 3.6 x 10^ 
K, then: 

Lx = 3.7 X iO-^^V^ / /n[4.7 x 10^° {T/u)]exp{-hu/kT)du ergs s"^ 

(Lang 1999), where V is the source volume and is the electron density. Taking the limits 
of integration to be over the soft X-ray band (0.2 to 2.5 keV) and Tc = 5 x 10^ K, (Model 
5), the above becomes: 

= 2 X 10^6 ( J^) ( ^^)' ergs s'^ 

^kpc^^ ^O.Olcm-3^ ^ 

A mean coronal density (between and Zstaii) of 0.035 cm~^ (Model 5, Table 3) in a 
rectangular region roughly the size of the modeled region of the plume (11 kpc x 13^ kpc) 
provides an estimate of the soft X-ray luminosity of — 4.6 x lO^'^ erg s"^. Using instead 
the density at Zstaii of 2.8 x 10~^ cm~^ gives Lj. — 2.9 x 10^^ erg s~^. The X-ray luminosities 
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of the hot gas in NGC 253 and M 82 are 2 x 10^^ erg s"^ and 2 x 10"^° erg s-^ respectively 
(Pietsch et al. 2000). Thus, the predicted X-ray halo for NGC 2613 should be observable 
below height, Zgtaih though is unlikely to have the smooth exponential distribution envisioned 
in this simple model. Note that NGC 2613 is more than an order of magnitude more massive 
than M 82 and therefore is more likely to retain any hot gas which might be around it. The 
values computed here will be lower, of course, if some other effect is contributing to the 
outflow. 



4.6. Blowout of Disk Gas 

The results of Table 3 show that it is feasible to transfer large masses to high galactic 
latitudes via buoyancy in a postulated X-ray corona. This drastically reduces the computed 
input energy requirements since it is no longer necessary to eject large masses to high alti- 
tudes. It is only necessary to achieve blow-out through the thin HI disk of NGC 2613. The 
conditions required for blow-out have been investigated by a variety of authors (see Tomisaka 
1998, for example) but energy requirements are typically of order 10^^ ergs, rather than the 
1q56 gj.gg ^-^^^ would normally be required for the -307 km s~^ feature. 

While the details of the interaction between the HI and hot outflowing gas are beyond 
the scope of this paper, it is important to ask whether there would originally have been 
sufficient HI in the disk from which the 8 x 10^ M© in the -307 km s~^ feature could have 
been swept up. Assuming cylinderical geometry, a disk region of radius, 6.5 kpc (Sect. 3.4) 
and using the modeled density and thickness of the disk (Table 2), the available HI mass is 
~ 3 X 10^ Mq suggesting that ~ 26% of the disk mass is swept up. This fraction would be 
higher for cone-like outflow. 

We should also consider whether the mass outflow rate is consistent with that of a 
galactic fountain. The estimated HI mass outflow rate (Sect. 3.4, Sect. 4.3) is ~ 1 to 2 Mq 
yr^^. Thus, the combined HI + hot gas outflow rate will be of order several Mq yr^^. Collins 
et al. (2002) have estimated the mass flow rate for the diffuse ionized gas (DIG) component 
in a galactic fountain, assuming a ballistic model of gas clouds. They find global values of 
M = 22 v^/y/0.2 Mq yr-i for NGC 891 and M = 13 ^/jv/^ M© yr'^ for NGC 5775, 
where fy is the filling factor. If these global values apply to a 20 kpc radius disk, taking 
fv — 0.2 and scahng to the disk area of the -370 km s~^ feature yields mass outfiow rates 
of 1.4 and 2.3 Mq yr^^ which are comparable to what we estimate, above. It is not yet 
clear, however, whether these galactic fountain values can be directly scaled to the relevant 
regions of NGC 2613. It may be that some additional source of pressure is still required, for 
example, magnetic fields in the form of a Parker instability. If such fields continue to rise 
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into the corona, the inclusion of this magnetic pressure would relax (i.e. lower) the density 
or temperature requirements internal to the plume (Table 3). 

5. Conclusions 

New VLA D array data of NGC 2613 have been combined with previous higher resolution 
observations (Chaves & Irwin 2001) to show a more extensive HI distribution than previously 
observed. The galaxy is now seen to have a tidal tail on its eastern side due to an interaction 
with its companion, ESQ 495-G017, to the north-west. 

The three-dimensional HI distribution in NGC 2613 has been modeled following Irwin 
& Seaquist (1991) and Irwin (1994), a method which allows the volume density distribution 
to be determined as well as the scale height and inclination to be disentangled. We find that 
the inclination of the galaxy (79°) is on the low end of the range given in Chaves & Irwin 
(2001) and the model now shows that there is no HI halo in NGC 2613. Rather, the global 
HI distribution is well fit by a thin disk of exponential scale height, = 188 pc. The use of 
such a model is very important in drawing conclusions about the presence or absence of a 
global halo in a galaxy of this inclination. Previous reports of a lagging halo from PV slices 
alone can largely be attributed to projection against the background disk. 

While there is no significant global HI halo in NGC 2613, there are more discrete disk- 
halo HI features than previously detected and these HI features achieve extremely high 
latitudes. Even though a tidal interaction is occurring, we suggest that most of the dis- 
crete kpc-scale features have been produced internally rather than from impacting clouds, 
although we do not rule out the possiblity of the companion galaxy having some indirect 
effect (e.g. triggering instabihties) . The presence of many discrete features may be related 
to the fact that the global HI disk is thin, favouring blow-out. The observed z heights are 
quite remarkable (e.g. up to 28 kpc). The -307 km s^^ feature, in particular, below the plane 
on the advancing side, reaching 22 kpc in z height and of total mass, (8 ± 2) x 10'^ Mq, is 
very obvious and well-defined in PV space. Its center is likely close to its projected radius 
of 15.5 kpc and it extends over a large (± 7 kpc radius) projected galactocentric radius. If 
this feature has achieved its z height as a result of internal processes, then extremely large 
energies are required, ~ 10^^ ergs. 

Given the very high input energies required for the -307 km s~^ feature, its resemblance 
to smaller buoyant features (cf. the Galactic Mushroom, Enghsh et al. 2000), and the 
fact that X-ray halos are being found around an increasing number of star forming spiral 
galaxies, we have carried out a feasibility study as to whether this feature can be interpreted 
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as an adiabatic buoyant plume. The observed HI would be carried out by a hot, low density 
outflowing gas and, after having achieved blowout, would rise through a hot pre-existing 
X-ray corona. A reasonable example (Model 5), has a mean plume temperature and density 
of 1 X 10^ K and 5.5 x 10^^ cm~^ rising into a hot isothermal corona of temperature and 
mean density, 5 x 10® K and 0.035 cm"^, respectively. These conditions produce a stall 
height of 11 kpc which is where the observed plume widens in velocity and position space. 
The coronal density at the stall height is 2.8 x 10~^ cm~^. The maximum outflow velocity 
in this model is 290 km s~^ and it reaches the stall height in 5.4 x 10^ yrs. This model shows 
that, even with buoyancy alone (and there may be additional sources of pressure such as 
magnetic fields), HI can reach these extreme z heights. 

The advantage of such a model is that the energy requirements from the initial event 
are drastically reduced to being only what is required for blowout, a reduction of several 
orders of magnitude. The energy is largely extracted from the gravitational potential of the 
galaxy rather than the initial event within tlie disk. The behaviour of the plume should 
sample the parameters of the X-ray corona. The predicted X-ray luminosity suggests that 
the corona should be observable at heights below the stall height although we expect that 
the distribution of X-ray emission may not be as smooth as assumed by the model. 

JI wishes to thank the Natural Sciences and Engineering Research Council of Canada 
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Fig. 1. — VLA C array HI Column density map of NGC 2613 superimposed on a Digitized 
Sky Survey optical image, from Paper I. Extensions are marked Fl through F6. 
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Fig. 2. — (a) HI total intensity (zeroth moment) map of NGC 2613 from the combined C+D 
array data, rotated so that the major axis is horizontal The central positions of NGC 2613 
and its companion, ESO 495-G017, are marked with crosses. The map was formed by first 
smoothing the data by 3 channels in velocity and 11 pixels (~ 1.4 beams, 1 pixel = 5") 
spatially, applying a cut off of 0.35 mJy bearn"^ to the smoothed, lower rms cube, and 
then summing the unsmoothed data excluding those points that fell below the cutoff. The 
greyscale ranges from to 3000 Jy beam~^ m s~^ and contours are at 20, 30, 50, 90, 150, 
250, 500, 1200, 2150, 2800 and 3200 Jy beam"^ m s"^ (1 Jy beam-^ m s"^ = 7.375 x 10^^ 
cm~^). Vertical numbered lines indicate where the slices were taken for the PV plots of Fig. 
3 and the features, Fl to F6 are again labelled, (b) The intensity-weighted mean velocity 
field (first moment map). Contours, in units of km s~^ are marked. This map was made in 
the same fashion as (a) but the cutoff was set to 0.45 mJy beam"-*^. 
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Fig. 3. — Position- velocity plots corresponding to the numbered vertical slices shown in Fig. 
2. Solid curves represent the data and red dashed curves represent the best fit model (see 
Table 2). The velocity is with respect to the systemic velocity. Contours are at 0.65 (1.5a), 
0.85, 1.1, 1.5, 3, 6, 10, 15, and 22 mJy beam^^. Negative contours (not shown) arc randomly 
distributed over the images. The last two panels show an average over 250" for the receding 
and advancing sides, respectively. For these, contours are at 0.2, 0.35, 0.6, 0.9, 2, 3, 4.5, 6, 
10, and 15 mJy beam~^ and the rms noise is ~ 0.16 mJy beam~^. 
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Fig. 4. — Sequence of 4 channels on the advancing side of the galaxy showing the -307 km/s 
feature with contours and greyscale chosen to highhght this feature. Contours are at -0.9, 
0.9 (2(t), 1.2, 1.5, 1.9, 3, 6, and 10 mJy beam~^ and the greyscale ranges from 0.45 (la) to 
10 mJy beam~^. Velocities are marked at the top right of each panel. 
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Fig. 5. — Zeroth moment map of the residuals (data cube - best fit model cube) 1 for all 
points > 3o". Solid curves represent excess data and dashed curves represent excess model, 
with contours at -800, -600, -400, -200, -100, 100, 200, and 400 Jy beam"! m 8"^ 
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Fig. 6. — Velocity as a function of 2; height for the five buoyancy models, Model 1 (crosses), 
Model 2 (diamonds). Model 3 (circles). Model 4 (boxes) and Model 5 (points), whose pa- 
rameters are given in Table 3. 



-29- 



Fig. 7. — z height as a function of time for the five buoyancy models, Model 1 (crosses), Model 
2 (diamonds). Model 3 (circles). Model 4 (boxes) and Model 5 (points), whose parameters 
are given in Table 3. The integration of these curves have been halted at Zstaii- 
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Table 1. Observing and Map Parameters. 



Parameter D Array C+D Array 

No. velocity channels 63 63 

Velocity Resolution (km s~^) 20.84 20.84 

Total Bandwidth (MHz) 6.25 6.25 
Synthesized Beam 

major X minor axis (" X ") @ PA (°) 96.5 X 42.9 @ -21.6 47.1 x 32.1 ® -8.2 

Rms noise/channel (mjy beam~l) 0.93 0.45 

Rms noise/channel (K) 0.14 0.18 
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Table 2: Modeled Kinematic and Density Parameters for Combined C+D Array Data. 



Parameter 


Whole Galaxy" 


RA (J2000) (h m s) 


08 33 22.8 ± 0.3 


DEC (J2000) (° ' ") 


-22 58 29 tt 


PA (°) 


114.2 t'oi 


in 


79.2 


Ysys (km s"^) 


1663 tl 


Vrnax (km S"^) 


304 ± 4 


R-maa; ( ) 


112 +^ 

-'--'-^ -59 


m 


U.Oi _Q 35 


(Ty (km s~^) 


17 ± 5 


no (cm"^)'' 


43 +° °^ 

'-'■^"J -0.26 


Ro (") 


140 ± 5 


Do (") 


95 +\ 


(") 


81 tl' 


He (") 


1 5 +2-0 

-■-■"3 -0.1 



"Errors reflect the variation in the parameter that results when the C array and D array data are modeled 
separately, and so throughout. 

''The density is the mean over the physical region corresponding to a 5" by 5" by 20.84 km s~^ (RA, DEC, 
Velocity) voxel. 
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Table 3: Buoyancy Model Parameters. 



Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


T, (10^ K) 


2.5 


2.5 


2.5 


2.5 


5.0 




0.74 


0.74 


0.50 


0.50 


0.74 




1.4 


1.4 


2.1 


2.1 


2.8 




100 


20 


100 


20 


10 




0.15 


0.15 


0.15 


0.15 


0.15 


ni{z^) (cm-3) 


0.0015 


0.0075 


0.0015 


0.0075 


0.015 


Z.tall (kpc) 


11.9 


7.4 


17.5 


11.0 


11.1 


n^^ (cm-3) 


0.015 


0.025 


0.016 


0.026 


0.035 


ndzstaii)^ {l^-^ cm-3) 


0.031 


0.76 


0.036 


0.80 


2.8 


Wi'' (10-3 cm-3) 


0.27 


2.1 


0.28 


2.2 


5.5 


(10^ K) 


6.8 


2.0 


6.9 


2.0 


1.2 



Note. — See Sect. 4.3 for an explanation of the parameters, a) computed from Eqn. (3) and Eqn. (1) 
using = 188 pc (Table 2). h) The coronal density at z* (the top of the thin disk) is computed from the 
parameters of the HI distribution at the position of the plume (see Table 2). c) The mean coronal density 
is determined from Eqn. (2), over the region from z^, to Zstall- d) The coronal density at height, Zstall from 
Eqn. (2). e) The mean plume density is determined from Eqn. (7), over the region from 2;* to z^taU- f) The 
mean plume temperature is determined from Eqn. (6), over the region from z^, to Zstall- 
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